Photon emission in neutral current interactions with 
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Abstract. We report on our study of photon emission induced by iiy ~ 1 GeV (anti)neutrino neutral cuiTent interactions 
with nucleons and nuclei. This process is an important background for Vf appearance oscillation experiments. At the relevant 
energies, the reaction is dominated by the excitation of the A(1232) resonance but there are also non-resonant contributions 
that, close to threshold, are fully determined by the effective chiral Lagrangian of strong interactions. We have obtained 
differential and integrated cross section for the (anti)neutrino-nucleon scattering and compare them with previous results. 
Furthermore, we have extended the model to nuclear targets taking into account Fermi motion, Pauli blocking and the in- 
medium modifications of the A properties. This study is important in order to reduce systematic effects in neutrino oscillation 
experiments. 
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INTRODUCTION 

A good understanding of (anti)neutrino interactions with nuclear targets in the E^ ^ \ GeV region is vital to reduce 
the systematic uncertainties in oscillation experiments aiming at a precise determination of neutrino properties. One 
of the possible reaction channels is photon emission induced by neutral current interactions (NC7), which can occur 
on single nucleons 

v{v)N^v{v)yN, (1) 



and on nuclear targets 
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with incoherent [Eq. (2)] or coherent [Eq. (3)] reaction mechanisms. It is also possible that, after nucleon knockout, 
the residual excited nucleus decays emitting /rays. This mechanism has been recently investigated [1] and shall not 
be discussed here. 

Weak photon emission has a small cross section compared, for example with pion production (NC;r), the most 
important inelastic mechanism. Indeed, while NC;r takes place predominantly via a weak interaction followed by a 
strong decay, in the case of NC7 one has a much weaker electromagnetic vertex instead of the strong one. In spite of 
this, NC7 turns out to be one of the largest backgrounds in v^ -^ Vf( V^ -^ Ve) experiments when 7's are misidentified 
as e^ from charge-current quasi-elastic scattering of Vf.(Ve). 

This is precisely the situation in the MiniBooNE experiment, where the gamma background is estimated from the 
measured NC;r*' rate assuming that it comes form radiative decay of weakly produced resonances, mainly A.^f Ny. 
The experiment finds an excess of events with respect to the predicted background in both v and v modes. In the v 
mode, the data are found to be consistent with v^i — > Vg oscillations and have some overlap with the LSND result [2]. 
In contrast, the reconstructed energy distribution of the e-like events in the v mode is only marginally compatible 
with a simple two-neutrino oscillation model, exhibiting an unexplained excess of events for E^ < 475 MeV [2, 3]. 
While several exotic explanations for this excess have been proposed, it could be related to unknown systematics and 
backgrounds. In the kinematic region where this anomaly is observed, NC7 is the second largest background behind 
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FIGURE 1. Set of Feynman diagrams of the model 

NCtt". In view of this, it is impoitant to study this process in detail, using the well developed framework of hadronic 
and nuclear physics. 

The first step in this direction was performed in Ref [4], where reaction (1) was studied with a microscopic model 
developed in terms of hadronic degrees of freedom: nucleon, A(1232) resonance and mesons. With this model, the 
NC7 events at MiniBooNE were calculated to be twice as many as expected from the MiniBooNE in situ estimate. 
The conclusion was that NC7 events give a significant contribution to the low-energy excess of e-like events [5]. 
However, in Ref. [5] the nuclear target C^C) was treated as an ensemble of nucleons, neglecting all nuclear-medium 
corrections. Furthermore, an energy independent and rather high efficiency correction compared with the presently 
available figures [6] was assumed in the analysis. A contrasting result, much closer to the MiniBooNE estimate, was 
recently obtained in Ref [7], based on the chiral effective field theory of nuclei [8], phenomenologically extended to 
the intermediate energies {Ey ~ 1 GeV) of the v flux at MiniBooNE. Our approach, described in the next section, has 
several ingredients in common with these previous works but also quantitative differences. 

FORMALISM 

First, let us discuss NC7 on the nucleon. The amplitude can be cast as 

.^,. = ^e;^''^a{p')r^"u{p)la. (5) 

It is proportional to the Fermi constant Gp and the electric charge e; e*'''' denotes the photon polarization vector; la, 
the Standard Model neutral current for (anti)neutrinos, is contracted with the hadronic current mF'^^m, which is specific 
for the different reaction mechanisms. The later is determined by the set of Feynman diagrams shown in Fig. 1 : direct 
and crossed nucleon-pole and A-pole terms, and pion-pole term 

PM«^rr + rr+rr. (6) 

For the nucleon-pole diagrams, F^ takes the form 

where Dn{p) = {j/ —M)^^ is the nucleon propagator; qis the 4-momentum transferred to the nucleon and qy, the one 
of the outgoing photon. 

JEMi-ly) = tFm-^(y^'qyvF2{0), (8) 

where fi 2 and Fi 2 are the vector NC and EM form factors, respectively. For the axial form factor 

we take 

^^(eV..(i-|) ' (10) 



with gA = 1 -267, the axial coupHng, and M^ = 1 .016 GeV [9], neglecting the strange part of the axial form factor F^ , 
as well as in the vector NC ones. 

The most prominent contribution to the cross section arises from the weak excitation of the A(1232) resonance 
followed by its radiative decay. For direct and crossed A terms one has 

r^" = jll,{p'Ay)Di,{p + q)J"Nc{PA)+j'N''c{p', -q)DUp' ^VTm^P^ -^y) > d D 



where !"■" = 70 {J"^) 7o and the A propagator is given by 

D^5a{p) = 



Age 



p^-Ml + iMiJTi^ip^) 



(12) 



A^fj is the A^ — A projector and T/^lp^) the A (energy dependent) width dominated by the A — > A^ TT p-wave decay. The 
vertices J^^ ^^ can be written in the most general form as 
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yp-qyp'') 


75, 


(14) 



in terms of the EM vector and NC vector and axial A^ — A transition form factors Cf^ = —Cf, Cf = —{1 — 2sin ©iv)^'^ 
and Cf — —Cf, respectively; C™ can be related to the A^ — A( 1232) helicity amplitudes, for which, following Ref. [10], 
we adopt the parametrizations obtained by the global analysis of n photo- and electro-production data with the unitary 
isobar model MAID [ 11 ] . In the axial sector we assume C^ —Q and Cf — —C^ /4 for the subleading (in a q^ expansion) 
form factors, while for the dominant C5 we take 



Ci{q^)^C^{Q){\ 



q 



0.93 ± 0.07 GeV fixed in a fit to v^ <i 



(15) 



IJ.-A++n BNL and ANL data [12]. The 



with C^ (0) = 1 .0 ± 0. 1 1 and Ma 

error in Cj (0) is the main source of theoretical uncertainty in our predictions. 

The first four diagrams in Fig 1 account for the leading contributions close to threshold. The inclusion of form 
factors allow to extend the model to the higher energies of interest for neutrino experiments. The last, pion-pole, term 
is of higher order and should be small. It is indeed found to be negligible compared to the mechanisms discussed 
above. We assume that other higher order terms can be also neglected. 

The model for NC7 on the nucleon has been extended to the incoherent reaction (2) on nuclear targets. For this 
purpose we have adopted the relativistic local Fermi gas description, according to which the target nucleons have 

momenta up to a Fermi momentum defined locally Pp'^{r) — [3n^Pp,n{r)] as a function of the local density of 
protons and neutrons independently. The density distributions are based on empirical determinations in the case of 
protons and on realistic theoretical models in the case of neutrons. Final nucleons are not allowed to take occupied 
states (Pauli blocking). 

Furthermore, it is known that the properties of the A(1232) resonance get substantially modified in a nuclear 
environment. This nuclear effect can be taken into account by performing the following substitutions in the A(1232) 
propagator [Eq. (12)] 



Ma 
Ta 



Ma + RcEaCp), 
fA-2ImZA(p). 



(16) 
(17) 



The real part of the in-medium A selfenergy, Ea, receives an attractive (negative) contribution from the nuclear mean 
field, which is partially cancelled by an effective repulsive piece from iterated A-hole excitations. In view of this and 
for the sake of simplicity we take ReEA(p) ~ 0. The resonance decay width is reduced to Fa because the final nucleon 
in A ^ nN can be Pauli blocked but, on the other hand, it increases because of the presence of many body processes 
such as AA^ -^ NN, AN — > NNn and AA^A^ -^ NNN (collisional broadening). These new decay channels, which are 
accounted in ImEA, have been parametrized as a function of the local density in Ref. [13]. 



RESULTS 



12 



10 



vp->vpv 



Hills FF(v 
Hills FF(v: 




0.2 0.4 0.6 

Ev (GeV) 



12 



10 



vn->vnv 
0.89GeV<C5(0)<1.11 GeV 

V ^^m ^^m 





1 


1 




^^ 1 






1 


i 



0.2 0.4 0.6 0.8 

Ev (GeV) 



FIGURE 2. Integrated cross section for (v)vp^ (v)vp7 (left panel) and (v)vn^ (v)vny (right panel). 

We consider first NC/on single nucleons. The integrated cross sections as a function of the (anti)neutrino energies 
are given in Fig. 2. The A mechanism is dominant but, at Ey ^ 1.5 GeV, the cross section from nucleon-pole terms 
is only ^ 2.5 smaller than the A one. The eiTor bands are determined by the uncertainty in C^ (0) discussed above. 
As in other weak interaction processes, the different helicities of v and v are responsible for different interferences, 
resulting in smaller v cross sections with a more Unear energy dependence. The dashed and dotted curves are obtained 
with the assumptions of Ref. [4]. The small differences, mainly in the faster rise of the cross section can be explained 
by: a larger Cj (0) = 1 .2 vi 1 ± 0. 11 of the present work, a constant Fa = 120 MeV of Ref, [4] vs the energy dependent 
width used in this work, and an Ma = 1 .2 GeV vs 1 GeV taken here for the nucleon-pole terms. 

In Fig. 3 the cross section for reaction (2) on '^C is presented. The dashed lines are obtained by summing the 
contributions of 6 protons and 6 nucleons at rest, i.e. neglecting all nuclear effects. By taking into account Fermi 
motion and Pauli blocking, the cross section already goes down by more than 10 %. With the full model, including 
the modification of the A resonance, the reduction is ~ 30 %. A similar net effect is obtained in Ref. [7], although the 
reduction quoted for the direct A mechanism (50 %) is substantially larger than ours (^ 30 %). 
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FIGURE 3. Integrated cross section for incoherent NCy on a '^C target. The error band corresponds to the uncertainty in 



C^(0) = liO.Il. 



Finally, using the available information about the total number of protons on target (6.46 x 10 in v mode and 
11.27 X 10^" in V mode), the target mass (806 tons) and composition (CH2) of the MiniBooNE detector [2], as 
well as the (anti)neutrino flux determination [14] and the energy-dependent efficiencies for 7 detection [6] we have 
calculated the number of e-like events from photons according to our model (further details will be given elsewhere). 
The comparison to the MiniBooNE in situ estimate is shown in Fig. 4 as a function of the (anti)neutrino energy 
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FIGURE 4. Distribution of e-like events at MiniBooNE as a function of the reconstructed (anti)neutrino energy. The solid line 
is the prediction of the present model while the dashed line is the MiniBooNE estimate based on NCtt measurements. The error 
band corresponds to the uncertainty in Cj (0) = 1 ± 0. 1 1 . 



reconstructed assuming charged-current quasielastic scattering on a bound nucleon at rest 

^QE _ 2{M- EB)Ey - {El - 2MEb) 
" 2[{M - Eb) - Ey{l - cosBy)]' 



(18) 



The Ve and Vg components of the flux can be neglected as expected but the wrong sign (V|, in v mode and viceversa) 
neutrinos yield a sizable contribution to the spectrum, in particular v^ in v mode. We find that our results are in good 
agreement with the MiniBooNE determination. In spite of the quantitative differences in the models, we arrive at the 
same conclusion as in Ref. [7], namely that NCy cannot explain the observed excess of e-Uke events at low E^ . 
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